The Atomic Slide Puzzle: Self-Diffusion of an Impure Atom. 
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In a series of recent papera3~H van Gastel et al have presented first experimental evidence, based 
on a series of STM images, that impure. Indium atoms, embedded into the first, close-packed layer 
of a Cu(OOl) surface, are not localized but make concerted, long excursions. Such excursions occur 
due to continuous reshuffling of the surface following the position exchanges of both impure and host 
Cu atoms with the naturally occuring surface vacancies. Van Gastel et al have also formulated an 
original lattice-gas type model with asymmetric exchange probabilities, whose numerical solution is 
in a good agreement with the experimental data. In this paper we propose an exact lattice solution 
of several versions of this model. 
PACS No: 68.35.Fx, 05.40.Fb, 66.30.Lw, 07.79. Cz 



Surface mobility is usually believed to be localized in 
the vicinity of the terrace steps and kinks of crystal sur- 
faces. These sites form the natural locations for atoms 
to attach to or to detach from the terraces; the diffu- 
sion of adatoms along or between the steps is known to 
induce their roughening and is responsible for the two- 
dimensional diffusive motion of adatom and vacancy is- 
lands By contrast, the atoms within the close- 
packed terraces that are not in the immediate vicinity of 
steps have often been considered as immobile since they 
are tightly packed by their neighbors. Although some 
suggestions have been made that the surface vacancies 
may diffusive themselves within the surface layers (see, 
e.g. Ref. Q), up to a very recent time there were no 
experimental techniques allowing for the direct observa- 
tion of the vacancies' or atoms' diffusion within the close- 
packed surfaces. 

In a recent Letter and two accompanying papers van 
Gastel et al [^[jsj presented an indirect experimental ev- 
idence that, remarkably, the atoms of the close-packed 
surfaces do undergo themselves continuous random mo- 
tion. In a set of ingenious STM experiments, van Gastel 
et al [I D have managed to introduce impure (Indium) 
atoms within the first layer of a Cu(OOl) surface and were 
able to follow their positions by analyzing the series of 
consequtive STM images. They found that the embed- 
ded In atoms make concerted, long excursions within the 
first surface layer, which motion can be only explained 
(see Refs. [| and @,| for more details) by continuous 
reshuffling of the surface in a way which resembles a slide 
puzzle due to position exchanges of both the Cu and In 
atoms with the naturally occuring surface vacancies Q . 

Van Gastel et al also proposed a simple model 
describing the In atom dynamics. In this model, the au- 
thors considered a terrace on a stepped Cu surface as a 
finite square lattice all L x L sites of which except two 
are filled with Cu atoms; the In atom is initially placed 
at the lattice origin and the vacancy - at one of the adja- 
cent sites. Both In and Cu atoms move randomly along 



the lattice by exchaning their positions with the vacancy, 
whose random walk terminates (and may reappear again) 
as soon as it reaches the lattice boundary. 

A salient feature of the physical situation is that, due 
to the difference of the Cu — Cu and the Cu — In interac- 
tions, the vacancy being at the adjacent to the In atom 
site has a preference to exchange its position with the 
In atom, compared to three adjacent Cu atoms. The 
chemical specificity of atoms is then taken into account 
by introducing unequal hopping probabilities; that is, in 
case when one of atoms adjacent to the vacancy is the In 
atom, the vacancy exchanges its position with the latter 
or with one of three Cu atoms with the probabilities 
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respectively. For the chemical species under study and 
at the room temperature the asymmetry parameter e as- 
sumes the value e* = 3 — I0~^, according to the Embed- 
ded Atom Model calculations In case when all four 
adjacent to the vacancy atoms are the Cu atoms, the 
vacancy exchanges its position with any of them with 
the probability = 1/4. As a matter of fact, numerical 
simulations joj show that the difference in the exchange 
probabilities has a significant impact on the In atom dis- 
placements: the average length of excursions is 2.2 times 
larger than in the e = case. Monte Carlo simulations 
of this model produce the results which are in a good 
agreement with the experimentally obtained jump length 
distribution and, hence, confirm the vacancy-assisted 
mechanism of the In atom dynamics proposed in Ref. [|| . 

Exact solution of this model is known only in the limit 
e = (see also Ref. |0 for the analysis in case of ar- 
bitrary vacancy concentrations p). Van Gastel et al 
have, however, furnished an approximate, continuous- 
space solution for the jump length distribution function 
supposing that the boundary conditions associated with 
the vacancy creation/annihilation appearing at the ter- 
race edges can be modeled by introducing some, a priori 
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given vacancy life-time. Despite the fact that it is pos- 
sible to find such values of this heuristic parameter at 
which the obtained results for the jump length probabil- 
ity distribution agree quite well with the experimental 
ones, it is, however, difficult to claim that such a com- 
parison can provide meaningful values of other physical 
parameters. Clearly, an assumption that the vacancy can 
spontaneously disappear at any lattice point, not neces- 
sarily at the terrace edges, might be a useful simplifi- 
cation but is not controllable. On the other hand, it is 
well-known that in two-dimensions the continuous-space 
description may incur significant errors in the numerical 
factors, which will result in modified values of the phys- 
ical parameters. Lastly, the continuous-space approach 
is hardly appropriate for studying the probability dis- 
tributions on scales of order of several lattice spacings 
only. Consequently, for meaningful interpretation of the 
experimental data an exact lattice solution of the model 
is highly desirable. 

In this paper we present an exact lattice solution of sev- 
eral versions of the model introduced by van Gastel et al 
1^-^ for arbitrary values of the parameter e, — 1 < e < 3, 
where positive (negative) values of e correspond to the 
physical situation in which the In atom has a preference 
(reduced probability) for the position exchanges with the 
vacancy. More specifically, we consider a situation with 
a single vacancy performing the just described random 
walk on a finite lattice with periodic boundary condi- 
tions, which mimic annihilation/creation of the vacancy 
at the terrace boundaries and also solve the case with a 
single vacancy and L = oo, which is appropriate for eval- 
uation of the transient dynamics. We note that in view 
of a very low concentration p of the vacancies this 
transient regime may be very important and may persist 
over a wide time range. In both situations we determine 
exactly the leading long-time asymptotical behavior of 
the mean-square displacement (MSD) and the diffusivity 
of a single In atom, as well as of the probability that the 
In atom appears at position X at time moment t. All our 
results will be presented in dimensionless form as func- 
tions of the discrete time n and for the lattice spacing set 
equal to unity. Dimension-dependent forms can be ob- 
tained by rescaling n t/At and X — X/tr, where At 
(« 10~®s) is the typical time interval between successive 
jumps and a 2.5 A) is the intersite distance Q. 

We begin with more precise definition of the model. 
Consider a two-dimensional, periodic in both xi and X2 
directions, square lattice of unit spacing all L x L sites of 
which except two (a vacancy and an In atom) are filled 
by Cu atoms (see Fig.l). The In atom is placed initially 
at the origin, while the vacancy is at position with the 
vector Y = e_i. Here, e^, i/ G {±1,±2}, denotes unit 
lattice vectors; ei (e_i) is the unit vector in the positive 
(negative) xi-direction, while 62 (e_2) is the unit vector 
in the positive (negative) a;2-direction. 

Next, at each tick of the clock, n = 1,2,3,..., a va- 
cancy exchanges its position with one of four neighboring 
atoms according to the probabilities defined in Eqs.(f]). 



That is, it has a preference (or, on contrary, a reduced 
probability for e < 0) to exchange its position with the 
In atom, when this one is at the neighboring site, and 
no preferency for exchange with any of atoms when all 
four adjacent to the vacancy sites are occupied by the 
Cu atoms. Hence, the vacancy moves randomly displac- 
ing the atoms in its path, including the In atom; apart 
from four "defective" sites adjacent to the In atom, it 
performs conventional symmetric random walk. 




Fig. 1. An illustration of the vacancy-assisted diffu- 
sion mechanism of the In atom (black sphere) on a 
lattice filled with the Cu atoms (grey spheres). The 
arrows depict the asymmetric hopping probabilities of 
the In and Cu atoms. 



Let Pn(X) be the probability that the In atom, which 
starts its random walk at the origin, appears at the site 
X at time moment n, given that the vacancy is initially 
at the site e_i. Then, following the approach of Ref. 
pO| , we write down the equation obeyed by P„ (X) as: 

+ 00 +00 +OC +CXD 

p—l mi — 1 nip — l mp-|-i— i^i h>p 

1- E ^:,(0|-e,jj F„,(0|e,Je_i) X 

X^e.,+...+e.,,X + '5x,0 ^1 - E (0| e_i ) ^ , (2) 

where i<"„(0|e^) is the probability that the vacancy, which 
starts its random walk at the site e^, arrives at the origin 
for the first time at the n-th step, while F„(0|ey|e^) 
stands for the conditional probability that the vacancy, 
which starts its random walk at the site e^, appears at 
the origin for the first time at the n-th step, being at the 
[n — l)-th step at the site e^. 
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General solution of Eq. can be written down in form 
of the following contour integral (see, Refs. and [p^): 



and 



(3) 



where the contour of integration C encircles the origin 
counterclockwise, (k-X) denotes the scalar product, and 

(4) 

In Eq.(^) I?^(k) denotes the following determinant: 

I?e(k)^det(I-T5(k)), (5) 

the elements ( T^(k) ) of the 4x4 matrix T^(k) obey: 

(T^(k))^^ = exp(i(k.e,)) F4(0|e,|e_^), (6) 

the indices v and fi assume successively the values 
1,-1,2,— 2; F^{0\ei,\e^^) is the generating function of 
the conditional first-visit probability F„(0|ei,|e_^), i.e. 



F5(0|e,|e_^) = J2 K(0|e,|e_^) C, 



(7) 



71=0 



and the matrix U^^'' (k) is given by 

^_^^^.(k.e,j ^(i_e-'(''-'))(I-Te(k))-i. (8) 

In turn, the MSD of the In atom is defined as 

dC d^P^ik) 



X2 = 4Dnie)n = 



2ti 



(9) 



|k|=0 



where ^^(k) is determined by Eq.(|^). 

We note now that as far as we are interested only in 
the leading large-n behavior of P„(X) and of the MSD, 
we may constrain ourselves to the study of the asymp- 
totic behavior of the generating function p(tr) ^j^. ^-j ^j^g 
vicinity of its singular point nearest to ^ = 0. Omitting 
the details of straightforward, but rather tedious calcu- 
lations, we find that the leading, in the combined k — > 
and <^ ^ 1~ limit, behavior of ^'^(k) follows 



P.(k)=^i(e)k'+^2(e)(l-0 + ---, 



(10) 



where 



(^-2)(^ + 2£2)(6-3)2 ^ 

J'^\€) — — ■ X 

^ ^ 2(57re- 16e - Stt) 

(-14eL2 -f QntL^ - 77re + QiiL^ - 6^^ 



Btt) 



((8e - 37re - 37r)L2 + 47re)2 



(11) 



^ . . 2(7r - 2)(3 - £)((3 - 6)L2 + 5e - 3) 
3(l + e)(16e-57re-h37r) 
((14e - 67re - Gtt -h 6)1^ + 7^ + 37r)2 
^ ((Be - 37re - 37r)L2 + At:^)^ 



(12) 




Fig. 2. Plot of D{e), Eq.(^), versus the parameter e 
for different values of the terrace width. The solid line 
corresponds to L = 15, the dashed - to L = 40 and 
the dotted - to L = 600, respectively. For notational 
convenience D{e) is multiplied by the terrace area L^. 



Then, substituting these expressions into Eqs.(P and 
(^), and noticing that 



^C/f)(k)P5(0|eJe_i)=-^i(e) k^ 



(13) 



we find that as n —^ oo, the diffusivity Dn{e) of the In 
atom, Eq.(0), tends to a constant value 



Die) 



3(l + e)(3-e) 



t((3- 



e)L2 



be 



(tt + 2L2) 



(^(67re - 14e + 6tt - 6)L^ - Tire 



37r 



(14) 



Note that D{e) appears to be a non-trivial, non- 
monotoneous function of the parameter e (see. Fig. 2); 
D{e) attains a maximal value when e = ec ~ 3 — 
12V7r — 2/L, and is exactly equal to zero for e — — 1 
and e = 3. In the case e = — 1, the In atom does not 
move at all, since it is not allowed to change its position 
with the vacancy, while in the e = 3 case it gets localized 
because at each time step it is forced to exchange its po- 
sition with the single available vacancy. Note that such 
a "localization effect" for e = 3 is, of course, specific to 
the situation with a single available vacancy; as a matter 
of fact, the indium atom can get delocalized already if 
a second vacancy is present. Consequently, the result in 
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Eq.(p^ concerns only the leading at p = behavior. For 
finite but small p, there should be a correction term to 
Eq. ( p^ proportional to . On the other hand, in view of 
a very small concentration of the naturally occuring va- 
cancies , observation of such correction terms (non-zero 
values of D{e = 3)) would require very large observation 
times. Note also that the ratio D{e = e*)/D{e = 0) is 
approximately equal to 1.8 for terrace width L = 400, 
which shows that, indeed, asymmetric exchange proba- 
bilities have also a strong effect on the In atom diflfusivity, 
not only on the average excursion length [||. 

Lastly, following Ref. we get from Eqs.(||) and (||), 
and the asymptotical expansions in Eqs.(p^) to (p^, that 
for sufficiently large n, such that ln(i) <C n <C L^, the 
scaling variables = xi/D{e)n and ^2 — X2lD{e)n are 
distributed according to a Gaussian law: 



converges to a limiting, non-Gaussian form: 



P(a,C2) = (2^)-i exp (- (e? +^|)/2). 



(15) 



We turn now to the case L = 00, which will allow us to 
obtain a meaningful transient behavior. Again, omitting 
the details of calculations, we find that the leading, in the 
combined k — > and ^ ^ 1^ limit, behavior of 2?{(k) 
obeys: 



I?^(k)=.Fi(6) k2-.F2(e)(ln(l-0 



(16) 



where 

•^i(e) 
and 

•^2(e) 



^ (7r-2)(6-3)^(3^6-76 + 3^-3) 
(57re- 16e-37r)(37re-8e + 37r)2 ' ' 



87r(7r - 2)(e - 3)2(37re -Te + Sn- 3)^ 
3(1 + e)(57re - 16e - 37r)(37re - 8e + 3tt)^ ' 

(18) 



while the sum 



(k)F^(0|e^|e_ 



-1) still obeys 

Eq.(|13[) with J-"i(e) replaced by ^lie)- Consequently, we 
find that in the leading in n order the diffusivity of the 
indium atom obeys, for e < 3, 



Dnie) 



3(1 + e) 



ln(n) 



7r(37re — 7e -I- 37r — 3) 



(19) 



and Z?„(e = 3) = 0. Hence, Dn{e) is a monotoneously 
growing function of e for e < 3, and is discontinuous at 
e = 3. The ratio Dn{e — €*)/Dn{t — 0) is now of order of 
1.9, which shows that also for infinite lattices asymmetric 
exchange probabilities enhance the In diffusivity. 

Finally, similarly to the approach of Ref. |l^ , we derive 
from Eci s.(|3|) and (jj), and the asymptotical expansions in 
Eqs.(^6^ to (|l8|), the asymptotic form of the probability 
distribution i-',i(X). We find that as n ^ 00, the distri- 
bution function of two scaling variables 



m,2 



'47r[(37r- 7)e-t-3(7r- 1)]n^i/2 
. 3(1 + e) ln{n) 



2^1,2, 



(20) 



(21) 



where Kg is the modified Bessel function of the zcroth 
order. Note that for e = 0, the results in Eqs.(p^), ( [20| ) 
and (pl|), as well as Dn{e) in Eqs.(lJ) and (|ig|), coincide 
with the ones obtained in Ref. pC[ |. 

In conclusion, we have presented here an exact lattice 
solution of several versions of the model originally devised 
by van Gastel et al H |^ to describe dynamics of an im- 
pure. Indium atom within the first layer of a Cu(OOl) 
surface. We have evaluated, for arbitrary values of the 
asymmetry parameter e, the long-time asymptotical be- 
havior of the impure atom MSD, as well as the limiting 
scaling forms of the probability distribution. Our ana- 
lytical results can be useful for further interpretation of 
the experimental data on dynamics of impure atoms in 
close-packed surface layers. 
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